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We present a search for the standard model Higgs boson in final states with an electron or muon 
and a hadronically decaying tau lepton in association with zero, one, or two or more jets using data 
corresponding to an integrated luminosity of up to 7.3 fb~^ collected with the DO detector at the 
Fermilab Tevatron collider. The analysis is sensitive to Higgs boson production via gluon gluon 
fusion, associated vector boson production, and vector boson fusion, and to Higgs boson decays to 
TT, WW, ZZ and bb pairs. Observed (expected) limits are set on the ratio of 95% C.L. upper Umits 
on the cross section times branching ratio, relative to those predicted by the Standard Model, of 22 
(14) at a Higgs boson mass of 115 GeV and 6.8 (7.7) at 165 GeV. 

PACS numbers: 13.85.Rm, 14.80Bn 



I. INTRODUCTION 

The standard model (SM) of particle physics postu- 
lates a complex Higgs doublet field as the source of 
electroweak symmetry breaking, giving rise to non-zero 
masses of the vector bosons and fundamental fermions. 
The mass of the SM spin-zero Higgs boson, H, that sur- 
vives after the symmetry breaking is not predicted, but 
is constrained by direct searches at the LEP [l|, Teva- 
tron [^ and LHC ^ colliders, to be in the range 115 
- 127 GeV at the 95% C.L. Precision measurements of 
W and Z boson and top quark properties [4| indicate a 
SM Higgs boson mass, rriH = QGlji \y\- Over the Higgs 
boson mass range 115 < mn < 150 the branching frac- 
tions vary considerably, with H -^ bb {H — ^ t^t^) being 
the dominant (subdominant) decays for mH < 135 GeV 
and H -^ W^W~ {H -^ ZZ) becoming important for 
rriH > 135 GeV. Previous analyses by the DO and CDF 
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Collaborations have mainly focused on the decay modes 
H ^ bb in the low mass region and H -^ WW with both 
W bosons decaying to an electron or muon in the high 
mass region. 

A previous DO publication Q reported a Higgs boson 
search in the tau lepton pair plus two jets final state, with 
one tau decaying to a muon and the other to hadrons, 
using 1.0 fb~^ of data. The CDF collaboration has re- 
cently reported a similar search in the tau lepton pair 
plus at least one jet Q. In this Letter, we report the 
results of three searches involving the production of tau 
leptons that extend the previous results by adding more 
data, increasing the trigger efficiency, adding new search 
channels, and considering additional signal contributions. 
The final states used are: (i) ^t plus zero or one jet (de- 
noted /j.tO), {ii) ^T plus two or more jets (/j.t2), and (m) 
er plus two or more jets (eT2). The /irO, ^t2, and 6x2 
analyses use data collected with the DO detector [8] cor- 
responding to integrated luminosities of 7.3, 6.2 and 4.3 
fb^^ respectively. The erO final state is not considered 
here as it suffers from large background and brings little 
increase in sensitivity. 

The Higgs boson production processes considered are 
(i) gluon gluon fusion (GGF), gg ^ H {+ jets); {ii) vec- 
tor boson fusion (VBF), qq -^ qqH; (iii) associated vec- 
tor boson and Higgs boson production (VH), qq — > VH, 
where T^ is a VF or Z boson, and V -^ qq (or Z ^^ vv 



in the case of firO); and (iv) associated Higgs boson and 
Z boson production (HZ), qq — >• HZ, with H ^ bb and 
Z -^ TT. The GGF, VBF, and VH processes are further 
subdivided according to the Higgs boson decay, H — ?> tt, 
H — ^ WW, or (for the ^rO analysis) H —>■ ZZ, and 
these subchannels are denoted as GGF,-,-, GGFww or 
GGF zz, etc. The fractional decompositions of signal 
contributions expected from Monte Carlo (MC) simula- 
tions are shown in Fig. [1] for the Higgs boson production 
cross section and branching ratios, and the event selec- 
tion requirements, discussed below. 
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FIG. 1: (color online) Fractional yields for H signals from 
MC simulations as a function of Mh for (a) the /^rO and (b) 
£t2 analyses. The yields for each signal process are plotted 
as solid lines for H ^ tt decays and as dashed lines for the 
H ^VV decays. 



Tau leptons can occur either through direct decays of 
the Higgs boson (at low mass) or indirectly from H — > 
VV with V decays to ts (at high mass). The leptons 
may arise from r decay or (at high mass) directly from 
V decay. The (.t channel is more uniformly sensitive to 
Higgs boson production over the full allowed mass range 
than are the dedicated H ^ bb oi H ^f WW -^ €tvv 
analyses, thus improving the sensitivity of a combina- 
tion of searches, particularly in the intermediate mass 
region around 135 GeV. In the following, "r" represents 
a hadronically decaying tau and "lepton (^)" denotes e 
or /i. 



II. TRIGGER 

The /irO and /ir2 data were collected from the full 
suite of DO triggers. The main contributors were the in- 
clusive high transverse momentum muon, \x -\- jets, and 
fi + T triggers. The trigger efficiency is determined in a 
two-step procedure starting from the measurement of the 
efficiency for inclusive muon triggers. This is measured 
using Z — )■ fifi candidates and parameterized as a func- 
tion of muon transverse momentum (pr), pseudorapidity 
(?]), azimuthal angle (0), and instantaneous luminosity. 
We then determine the ratio of the yields of the full trig- 
ger suite relative to those for the inclusive muon triggers. 
For the firO analysis the ratio is parametrized as a func- 
tion of p^ while for the /zr2 analysis it is a constant. The 
efficiency for Z ^- tt events for the full suite of triggers 
varies between 80% and 95%, and is about 40% larger 
than for the inclusive muon trigger. 

For the 6x2 analysis, a set of calorimeter-based inclu- 
sive electromagnetic object triggers was used. The ef- 
ficiency of these triggers, obtained from an analysis of 
Z — > ee events selected with just one identified electron, 
is found to be about 85%. 



III. BACKGROUND AND SIGNAL SAMPLES 

The major backgrounds for the Higgs boson search are 
Z + jets, W + jets, ti, and multijet production (MJ) with 
misidentification of leptons or taus. Smaller backgrounds 
arise from boson (W, Z or 7) pair production and sin- 
gle top quark production. All but the MJ background 
are simulated using MC event generator programs and 
normalized to the highest available next-to-leading order 
(NLO) or next-to-NLO (NNLO) theoretical calculations. 
The MC simulations use the CTEQ6L1 parton distribu- 
tion functions (PDF) [9]. 

The Z + jets and W + jets event samples are generated 
by ALPGEN 10], interfaced to PYTHIA ll| which provides 
initial and final state radiation and hadronization of the 
produced partons. The pj. distribution is reweighted to 
agree with the DO measurement [I4I. The p^ is also 
reweighted for the ^t2 analyses using the reweighting fac- 
tors derived for the p^ distribution, multiplied by the 
ratio of the Pjf to the pf distributions as predicted in 
NNLO QCD [il|. For the £t2 analyses, the absolute nor- 
malization for the Z + jets and W + jets cross sections 
are taken from Ref. ^ using the MRST2004 NNLO 
PDFs 



15j . The same Z -I- jets normalization is used for 
the /xrO analysis but the W + jets normalization is de- 
rived from data as discussed below. 

We simulate ti and single top quark events using the 
ALPGEN and COMPHEP ilf)] generators respectively, with 
PYTHIA used to simulate hadronization effects. The nor- 
malizations are based on the approximate NNLO calcula- 
tions [l3l- The diboson events are generated by PYTHIA. 

Higgs boson production is simulated using PYTHIA, 
with normalizations taken from Ref. fig. We use HDE- 



CAY [1^ and TAUOLA \2d\ to obtain the branching frac- 
tions of the Higgs boson and tau lepton respectively. 

All MC signal and background events are input to a 
GEANT3-based f21'| simulation of the detector response 
and processed with the same reconstruction programs as 
used for data. Data events collected from random beam 
crossings are superimposed on the MC events to account 
for detector noise and pileup from additional pp collisions 
in the same or previous bunch crossings. Correction fac- 
tors are applied to the simulated events to account for 
the trigger efficiencies and for the differences between 
MC and data for the lepton, tau, and jet identifications, 
and for the energy scale and resolution of jets. 



IV. EVENT SELECTION CRITERIA 

Muons selected for this analysis are required to have 
hits in the muon chambers before and after the toroidal 
magnets and to be matched to a track in the tracking 
system with px > 15 GeV and |ry| < 1.6. Muon candi- 
dates are required to be isolated in both the calorimeter 
and the tracking system using the calorimeter transverse 
energy, Ei^°, in the annular cone 0.1 < TZ < 0.4 around 
the muon, where TZ = •^/(A?])^ -I- (A^)^, and the track 
transverse momentum sum, plf° — Ep^^'^'^, within a cone 
TZ < 0.5, excluding the px of candidate muon. For the 
/irO analysis, E'f° and plfP must be less than 15% of p^. 
For the fiT2 analysis, E!^° and Px° must be less than 2.5 
GeV. Muon candidates due to cosmic rays are rejected if 
the scintillation counters surrounding the detector indi- 
cate a time of arrival different by more than 10 ns from 
that expected for collision products. 

Electrons are identified using a likelihood variable, Ce, 
that uses as inputs the quality of the matching of the elec- 
tromagnetic (EM) shower centroid to a track, the fraction 
of energy deposited in the EM section of the calorime- 
ter (EMF), a measure of the probability that the energy 
deposit pattern in the calorimeter conforms to that ex- 
pected for an electron, E^°, and the separation along the 
beam axis of the electron track and the primary vertex 
(PV) [22|. The signal sample electrons are required to 
have £e > 0.85. Electron candidate tracks are required 
to have pt > 15 GeV and \r]\ < 1.1 or 1.5 < \r]\ < 2.5, and 
to impinge upon a module of the central EM calorimeter 
within the central 80% of its azimuthal range. 

The selection of hadronically decaying tau leptons is 
done separately for three types based on the number 
of tracks within a cone TZ < 0.3 and the number of 
EM subclusters found in the calorimeter using a near- 
est neighbor algorithm. Type-1, patterned on the de- 
cay r — >■ irvr, requires one track and no EM subclusters. 
Type-2, based on r — > p{tt TT^)i'r, requires one track and 
at least one EM subcluster. Type-3, motivated by the 
T -> 7r^7r^7rT(7^°)^'^- decay, requires at least two tracks 
with or without EM subclusters. We reject type-3 can- 
didates with exactly two tracks of opposite signs since 
their charge sign is ambiguous. The r transverse energy. 



E^ , is defined as the visible transverse momentum of the 
T decay products as measured by the calorimeter with 
appropriate energy scale corrections. The ratio of E^ to 
the sum of the tracks associated with the tau, p^^, is 
used to verify that the MC and data tau energy scales 
are the same. We require E^ > (12.5, 12.5, 15) GeV, 
p^f^ > (7,5,10) GeV, and {pt^/E:^) > (0.65,0.5,0.5) for 
T types (1, 2, 3). The leading (highest px) track for t ype - 
3 rs must exceed 7 GeV. A neural network, NN-r [23| . 
based on energy deposition patterns and isolation crite- 
ria in the calorimeter and tracking systems is constructed 
for each tau type to discriminate a r from a misidenti- 
fied jet. Lower bounds placed on NN,- at 0.9, 0.9 and 
0.95 for tau types 1, 2, and 3 select hadronically de- 
caying taus with good purity. For type-2 r leptons we 
discriminate taus from electrons using a second neural 
network, NN^/e, constructed using variables that char- 
acterize the longitudinal and transverse energy profiles 
in the calorimeter, the energy and position correlations 
between r tracks and calorimeter energy deposits, and 
isolation of the calorimeter energy. 

Jets are selected using an iterative midpoint cone al- 
gorithm [2J] with a cone size TZ — 0.5. We require at 
least two tracks associated with the jet that point to the 
PV. Jet energies are corrected to the particle level for 
out-of-cone showering, underlying event energy deposits 
and pileup, and the estimated missing energy in jets with 
identified scmileptonic decays of a hadron. The energy 
scale, resolution, and jet identification efficiency for MC 
jets are corrected to give agreement with data. For the 
quark-dominated MC samples {tt and diboson), there is 
an additional correction of the jet energy that accounts 
for the differences in the responses of quark jets and the 
dominantly gluon jets for which the jet energy scale cor- 
rection was obtained. The /xr2 and 6x2 analyses require 
at least two jets with |ryjot| < 3.4 and p!^^ > 20 (15) GeV 
for the leading (other) jet. The /zrO analysis imposes 
these jet px requirements as a veto to ensure that the 
selected samples have no events in common. 

The missing transverse energy, _^3^, is computed from 
the observed transverse energy deposits in the calorime- 
ter and is adjusted for the appropriate energy scale cor- 
rections for all objects and for isolated muons observed 
in the event. 

For the final event selection, all three analyses re- 
quire exactly one isolated lepton and a hadronic tau 
with opposite charges. The separations between all pairs 
of lepton, tau, and jet are required to be 7?. > 0.5. 
For the firO analysis, events are required to have only 
one T, and the smaller of the transverse masses, mx — 

2E^^ °"^t(1 — cos A0) (where "lepton" = r or /j, and 
A0 is the angle between the lepton andlf^x) must exceed 
25 GeV to suppress the Z-\- jets and MJ backgrounds, 
while retaining about 80% of the signal. For the eT2 anal- 
ysis, substantial backgrounds arise from Z -\- jets produc- 
tion with Z — >■ ee where an electron is misidentified as 
a type-2 r. To reduce these, we remove r candidates 



in the region 1.1 < \r]\ < 1.5 where the calorimetry 
has impaired electron identification. Further Z+ jets 
rejection is obtained by requiring type-2 t candidates 
to have NN^/g > 0.95 to suppress electrons that resem- 
ble the track + EM cluster signature. This cut retains 
more than 80% of type-2 ts while rejecting about 90% of 
the electrons [2^. We reject type-2 t candidates which 
point near the edge in (j> of an EM module in the central 
calorimeter where the EM response is impaired. In addi- 
tion, type-3 T candidates with EMF > 0.95 are excluded. 
The MJ background in the er2 analysis is suppressed by 
requiring S > 1, where 5 is a measure of the significance 
for J^T to differ from zero |26| . 



V. BACKGROUNDS DERIVED FROM DATA 

The MJ background arising from misidentification of 
leptons or taus by the detector reconstruction algorithms 
is difficult to simulate, so for each analysis, the MJ back- 
ground is taken from data. The general method for all 
analyses is similar: we define a sample of MJ-enriched 
events, A^, from which residual backgrounds simulated 
by MC are subtracted, to provide the shapes of the MJ 
kinematic distributions. The number of MJ events in the 
signal sample is obtained by multiplying the MJ yield in 
a signal- like sample A/" by a scale factor p^, obtained from 
the A^ sample for each of the tau types, i. The pi factors 
provide the estimate for the differences in the MJ back- 
ground normalization between the Af and signal samples, 
based on the A4 sample, and are in all cases within 10% 
of unity. 

For the /irO channel, the sample AA is obtained by re- 
quiring mTip^ST) < 30 GeV and NN^ < 0.2, and the pi 
are the ratios of isolated to non-isolated lepton events in 
A4, and are parameterized as a function oi p^, Nj^ts, ^t, 
and pt^. These factors scale the MJ fraction of the sam- 
ple Af, selected as for the signal sample except that the 
muon is required to be non-isolated, to obtain the MJ 
normalization in the isolated lepton signal sample. An 
alternate MJ-enriched sample is defined by NN^. < 0.2 
and mT{p,IpT) < 30 GeV, in which the r and p have 
the same charge sign, for estimating the MJ background 
uncertainty. 

For the pt2 analysis, the MJ sample Ai is obtained by 
reversing at least one of the muon isolation requirements 
and requiring 0.3 < NN^ < 0.8. The MJ fraction of this 
sample is 94% before the MC-simulated background sub- 
traction. The Pi factors are the ratio of opposite charge 
sign (OS) and same charge sign (SS) p — t pairs in A1 
and are used to scale the MJ component of the sample Af 
selected as for the signal sample except that we require 
SS p and t. The pi show no significant dependence on 
the kinematic variables. 

For the er2 analysis, A4 is obtained by requiring the 
electron to satisfy an orthogonal loose electron selection, 
0<Ce< 0.85, and 0.3 < NN^ < 0.9. The MJ fraction of 
this sample is 96% before the MC-simulated background 



subtraction. The pi are obtained from the OS and SS 
A4 sample and are applied to the MJ component of the 
Af sample as in the /it2 analysis. The pi show no sig- 
nificant dependence on kinematic variables. Alternate 
MJ-cnriched samples, in which either the t or lepton se- 
lections (but not both) are reversed, are defined for es- 
timating the MJ background uncertainties in both 1x2+ 
analyses. 

For the prO analysis, the dominant background is from 
W + jets with the muon from W decay and a jet misiden- 
tified as a tau. Both the normalization of the W + jets 
sample and the misidentification probability are difficult 
to model adequately, so the simulation is corrected using 
a data-driven method [23 • The jet produced in asso- 
ciation with a W boson has a charge that is correlated 
differently with the W boson charge for quarks and glu- 
ons. Furthermore, the probability for a jet misidentified 
as a tau to have the same charge sign as its progeni- 
tor parton varies with NN,-. We determine a weight for 
W + jets MC events that depends on the charge corre- 
lation between the muon and recoil parton and on the 
value of NN^. 



VI. EVENT YIELDS 

The numbers of data and expected background events 
are given in Table U for the yurO, yur2, and 6x2 analyses. 



TABLE I: For each analysis channel, the number of back- 
ground events expected from MC simulated processes, MJ 
background, and observed data, for individual and sum of all 
tau types after preselection. "V+j" denotes W ot Z + jets 
and "DB" denotes diboson processes. 



r type 


tt 


W+i 


Zu+i 


Zrr+i 


DB 


MJ 


EBkd Data 








prO 


analysis 








type 1 


4 


234 


22 


11 


29 


39 


338 


340 


type 2 


19 


852 


94 


56 


108 


116 


1245 


1294 


type 3 


4 


678 


57 


19 


25 


67 


850 


839 


All 


27 


1764 


172 


86 


162 


223 


2433 


2473 








pr2 


analysis 








type 1 


13 


9 


5 


29 


2 


19 


76 


81 


type 2 


86 


57 


22 


159 


12 


59 


394 


418 


type 3 


13 


34 


4 


43 


3 


22 


119 


109 


All 


112 


100 


31 


231 


16 


99 


589 


608 








er2 


analysis 








type 1 


2 


2 





6 


1 


6 


18 


10 


type 2 


14 


21 


14 


30 


2 


25 


106 


98 


type 3 


7 


16 


2 


11 


1 


15 


52 


59 


All 


24 


40 


16 


46 


4 


47 


176 


167 



VII. MULTIVARIATE ANALYSIS 

The expected number of events for Higgs boson sig- 
nal processes is small in comparison to the backgrounds 
shown in Table HI For example, the expected signal yields 



at iriH — 165 GeV are 5.2, 1.7 and 0.3 events for the firO, 
Ht2 and eT2 analyses respectively. The corresponding 
yields at tuh = 115 GeV are 0.9, 1.6 and 0.4 events. We 
thus employ multivariate techniques that utilize both the 
magnitudes of the variables and the correlations among 
them to separate the signal from the backgrounds. We 
choose well-modeled variables that have the capability 
to distinguish between at least one signal and one back- 
ground as shown in Table HIl Figure [2] shows distributions 
for representative variables that offer significant discrim- 
ination of signal and background for each of the chan- 
nels. For calculating the tt invariant mass shown in 
Fig- mjc), the .^T is apportioned to the neutrinos from 
the two postulated tau leptons by decomposing the I^t 
vector into components associated with the observed lep- 
ton and hadronic tau 12811. 



signals against the sum of all backgrounds, with all r 
types combined for Higgs boson masses 105 < mn < 200 
GeV in 5 GeV steps. The BDT output is binned in 15 
bins spanning — 1 < BDT < 1 with a non- uniform bin- 
ning to assure sufficiently small statistical uncertainty in 
the predicted backgrounds within any bin. We smooth 
the effects of signal MC statistics by averaging BDT dis- 
tributions for mn with the neighboring distributions at 
{ttih — 5) GeV and {niH + 5) GeV with weights of 50%, 
25%, and 25% respectively. Figure [3] shows the NN// 
distribution for the ^rO analysis at niH = 165 GeV and 
the averaged BDT distributions for the £t2 analyses at 
rriH — 150 GeV, where the sensitivities are maximal. 



VIII. SYSTEMATIC UNCERTAINTIES 



TABLE II: List of variables used in the multivariate discrim- 
inants for the firO and £t2 analyses. The variable _^t is the 
missing transverse energy computed from the jets in the event. 



Variable 



UtO ^r2H 



Lepton pt 
Tau Pt 
Leading jet pr 

]^T 

fi charge xrjfj, 

rjr 

It IpT invariant mass 

TT invariant mass f28j 

Dijet invariant mass 

/i T If^T invariant mass 

(.V transverse mass, mj, 

TV transverse mass, mJj, 

Minimum of m^ and mj- 

E|pt! of all jets 

Scalar pt sum of i, t, ]^t, jets 

Magnitude of vector pr sum oi £,t, J^t, jets 

Minimum y^s necessary for final objects x 

Number of jets x 

AT?, between leading jets x 

Ar; between leading jets x 

Asymmetry between _^t and^T x 

A(f> between £ and r x 

A9 between £ and r x 

A(j) between £ and ]^t x 

A(f> between r and_^T x 

A(f) between _^T from calorimeter and tracks x 

Cosine of angle between £ and beam direction x 

Minimum Stf) between ]^t and a jet x 

Missing Et significance, S x 

NN^ X 



X X 

X X 

X X 

X X 

X 

X 

X 

X 
X 
X 

X 
X 
X 

X 
X 
X 



The ^rO analysis uses neural networks ^29] (NN/f) 
trained to discriminate between all backgrounds and all 
signals for 115 < mn < 200 GeV in 5 GeV increments. 
Type-2 t samples are trained separately, while the t 
types 1 and 3 are combined for training to increase statis- 
tics. The NNh distributions are binned in 21 equal sized 
bins for < NN// < 1.05. The ^r2 and er2 analy- 
ses use boosted decision trees (BDT) [23I trained for all 



A large number of systematic uncertainties have been 
considered, typically broken down separately for each 
analysis channel, tau type, background or signal process, 
or Higgs boson mass. The luminosity and trigger uncer- 
tainties arc obtained from separate analyses of DO data. 
The lepton, tau, and jet energy scale, resolution, and 
identification uncertainties are obtained from special con- 
trol samples. Uncertainties in the MC-simulated back- 
ground cross section normalizations and shapes are ob- 
tained using theoretical uncertainties, and the extent to 
which special data samples enriched in each background 
process agree with MC predictions. The MJ background 
uncertainties are determined by comparing the alternate 
MJ-enriched samples with the results obtained with the 
nominal choice. Signal cross section uncertainties are ob- 
tained from theoretical estimates and include the effect 
of PDF uncertainties. For each source, the impact on the 
final variable (NN// or BDT) distribution is assessed by 
changing the nominal values of a parameter by ±1 s.d. 
Some of the uncertainties affect only the normalization 
of the final variable distribution and some also modify its 
shape. 

Table IIIII summarizes the systematic uncertainties. 
Many entries comprise several subcategories. For exam- 
ple, the jet reconstruction uncertainty includes the effects 
of jet identification, confirmation that the tracks within 
the jet arise from the PV, jet resolution and jet energy 
scale. Moreover, these elements of the jet reconstruc- 
tion uncertainty are computed separately for different 
background processes and hypothesized Higgs mass val- 
ues in each analysis channel. The dominant systematic 
uncertainties are due to the ^-l-jets and MJ backgrounds, 
with significant contributions from jet reconstruction and 
modelling for the £t2 analyses. 



IX. CROSS SECTION LIMITS 

We observe no excess of events over that expected from 
backgrounds in Fig. [31 We therefore obtain upper limits 
on the Higgs boson cross section for each analysis from 
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FIG. 2: (color online) Comparison of data, expected backgrounds, and total signal with the indicated scaling factors for (a) 
invariant mass of the /i, t, I^t system for the /irO analysis (signal shown for mn = 165 GeV); (b) dijet invariant mass for the 
fj,T2 analysis (signal shown for mn ~ 150 GeV); and (c) invariant mass of the tt system for the er2 analysis (signal shown for 
niH = 150 GeV), where the I^t contribution is apportioned to the e and t as discussed in the text. The "MET" in the labels 
refers to IpT- The it, single top and diboson backgrounds are shown together as "other". 
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FIG. 3: (color online) (a) NN_fr distribution for the /xrO analysis at ttih ~ 165 GeV; (b) BDT distribution for the ^t2 analysis at 
niH = 150 GeV; and (c) BDT distribution for the er2 analysis at niH ~ 150 GeV. The it, single top and diboson backgrounds 
are shown together as "other" . 



the final multivariate outputs using the modified frequen- 
tist method [30], using a negative log likelihood ratio 
(LLR) for the background only and signal+background 
hypotheses as the test statistic. For the /xrO analysis, 
each tau type is input separately to the limit setting cal- 
culation for Higgs boson masses from 115 to 200 GeV in 5 
GeV steps. The £t2 calculation uses the BDTs summed 
over tau type for niH values from 105 to 200 GeV in 5 
GeV steps, averaged over neighboring mass bins as de- 
scribed above. 

The impact of systematic uncertainties on the limits 
is minimized by maximizing a likelihood function [31] 
in which these uncertainties are constrained to Gaus- 
sian priors. The value of the Higgs boson cross section 
is adjusted in each limit calculation until the value of 
CLs reaches 0.05, corresponding to the 95% C.L., vifhere 
CLg — CLg+b/CLb and CLg+b (CLb) are the probabil- 
ities for the negative LLR value observed in simulated 
signal-f background (background) pseudo-experiments to 
be less than that observed in our data. The limits ob- 
tained are summarized in Table IIVI 



We combine the information from the three channels 
by recomputing the LLR and limits for the three analyses 
together, now also including the limits from the previous 
independent /iT2 analysis using 1 fb~^ [6|. In this calcu- 
lation, the systematic uncertainties across the different 
analyses are appropriately correlated (e.g. the Z + jets 
normalization for all channels is the same). The fully 
combined LLR distributions and the 95% C.L. limits as 
a function of mn are shown in Figs. |4] and [5j The com- 
bined limits are also shown in Table IIVI 

In summary we have searched for the SM Higgs bo- 
son in final states involving an electron or muon and a 
hadronically decaying tau. We set 95% C.L. limits on 
the Higgs boson production cross section which are 21.8 
and 6.8 times those expected in the SM for Higgs boson 
masses of 115 and 165 GeV. 
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TABLE III: The range of systematic uncertainties (in percent) 
for categories of their source. Each category generaUy summa- 
rizes several individual sources separated by analysis channel, 
tau type, and/or physical process. Those with "Type" indi- 
cated as "N" affect only the normalization of the final variable 
distribution. Those indicated as "S" also affect the shape of 
the final variable distribution. 
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Source 


Type 


Uncertainty 


Luminosity 


N 


6.1 


Muon trigger 


N 


5-9 


Electron trigger 


N 


2 


Muon reconstruction 


N 


2-3 


Electron reconstruction 


N 


4 


Tau reconstruction 


N 


4-14 


Jet reconstruction 


S 


2-10 


Jet modeling 


S 


0-7 


MC simulated backgrounds 


N 


5-12 


MJ background 


S 


10-50 


Signal cross sections 


N 


5-40 



TABLE IV: The ratio of expected and observed 95% C.L. lim- 
its on the Higgs boson cross section to the SM values for each 
analysis channel and the combination of all channels including 
that oi^. 



rriH 


firO 


^t2 


er2 


Combined 




exp 


obs 


exp 


obs 


exp 


obs 


exp 


obs 


105 


- 


- 


17.7 


18.5 


33.3 


58.9 


12.6 


17.1 


110 


- 


- 


19.3 


20.8 


34.3 


55.7 


12.9 


17.7 


115 


84.2 


106.4 


20.3 


26.3 


37.5 


55.1 


14.3 


21.8 


120 


42.9 


31.1 


19.2 


23.3 


40.5 


59.4 


13.7 


15.6 


125 


34.2 


37.5 


17.3 


19.5 


42.3 


64.9 


12.8 


15.7 


130 


25.2 


32.4 


15.9 


20.6 


44.2 


72.5 


11.5 


17.9 


135 


20.3 


20.3 


17.5 


15.2 


47.2 


82.5 


11.3 


11.8 


140 


16.7 


20.0 


18.7 


13.2 


44.7 


68.1 


11.1 


10.1 


145 


13.8 


13.3 


18.3 


12.9 


43.5 


54.2 


11.3 


9.8 


150 


11.9 


12.8 


17.9 


13.6 


45.4 


54.1 


10.8 


9.5 


155 


9.8 


12.9 


18.2 


13.2 


42.3 


57.5 


9.2 


9.0 


160 


8.2 


7.6 


19.1 


11.1 


33.9 


74.9 


8.4 


7.6 


165 


8.1 


7.8 


21.7 


11.2 


32.8 


69.8 


7.7 


6.8 


170 


8.5 


9.4 


21.3 


12.7 


35.2 


64.5 


8.5 


7.4 


175 


9.5 


8.6 


22.7 


11.4 40.7 


73.7 


9.6 


8.0 
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13.5 


22.1 


14.6 


45.5 


84.6 


11.4 


11.0 
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13.5 


12.1 


25.7 


19.8 


53.7 


90.8 


12.2 


9.7 
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16.5 


17.2 


29.5 


19.1 


58.8 


101.8 


14.6 


12.3 
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18.5 


18.7 


30.1 


20.9 


67.3 


110.4 


16.1 


15.3 
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19.2 
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28.9 


26.9 


69.3 


114.4 


19.8 


29.9 
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